Observational Constraints on a Variable Dark Energy Model 
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We study the effect of a phenomenological parameterized quintessence model on low, intermediate 
and high redshift observations. At low and intermediate redshifts, we use the Gold sample of super- 
nova Type la (SNIa) data and recently observed size of baryonic acoustic peak from Sloan Digital 
Sky Survey (SDSS), to put constraint on the parameters of the quintessence model. At the high 
redshift, the same fitting procedure is done using WAMP data, comparing the location of acoustic 
peak with that obtain from the dark energy model. As a complementary analysis in a flat universe, 
we combine the results from the SNIa, CMB and SDSS. The best fit values for the model parame- 
ters are fi m = 0.27to'o2 (^ ne P resen t matter content) and wo = — 1.45lg'go (dark energy equation 
of state). Finally we calculate the age of universe in this model and compare it with the age of old 
stars and high redshift objects. 



PACS numbers: 05.10.-a ,05.10.Gg, 05.40.-a, 98. 
I. INTRODUCTION 

Observations of the apparent luminosity and redshift of 
type la supernovas (SNIa) provided us with the main evi- 
dence for the accelerating expansion of the Universe [1,2]. 
A combined analysis of SNIa and the Cosmic Microwave 
Background radiation (CMB) observations indicates that 
the dark energy fills about 2/3 of the total energy of the 
Universe and the remaining part is the dark matter and 
a few percent in the Baryonic form [3-5]. 

The " cosmological constant" which was introduced by 
Einstein to have a static universe, can be a possible solu- 
tion for the acceleration of the universe [6]. The cosmo- 
logical constant in Einstein field equation is a geometri- 
cal term, however, it can be regarded as a fluid with the 
equation of state of w = —1. There are two problems 
with this scenario, namely the two dark energy problems 
of the fine-tuning and the cosmic coincidence. Within 
the framework of quantum field theory, the vacuum ex- 
pectation value of the energy momentum tensor diverges 
as fc 4 . A cutoff at the Planck energy leads to a cosmologi- 
cal constant with 123 orders of magnitude larger than the 
observed value of 1CU 47 GeV 4 . The absence of a funda- 
mental mechanism which sets the cosmological constant 
to zero or to a very small value is the cosmological con- 
stant problem. The second problem, cosmic coincidence, 
states that since the energy densities of dark energy and 
dark matter scale so differently during the expansion of 
the Universe, why are they nearly equal today? 

One of the solutions to this problem is a model with de- 
caying cosmological constant from the Planckian area at 
the early Universe to a small enough at the present time. 
Dolgov (1983) proposed a massless non-minimally cou- 
pled scalar field to the gravity with a negative coupling 
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constant to solve this problem [7]. However, this model 
provides a time varying Gravitational constant which 
strongly contradicts the upper limits from the Viking 
radar range [8] and lunar laser ranging experiments [9] . A 
non-dissipative minimally coupled scalar field, so-called 
Quintessence model can also play the role of time varying 
cosmological constant [10-12]. The ratio of energy den- 
sity of this field to the matter density increases by the ex- 
pansion of the universe and after a while the dark energy 
becomes the dominated term in the energy-momentum 
tensor. Tuning the parameters of this model it can pro- 
duce the value of A-term both for the early universe and 
the present time. 

One of the features of the quintessence model is the 
evolution of the equation of state of dark energy dur- 
ing the expansion of the Universe. Various models de- 
pending on the potential for the scalar field as k-essence 
[13], Tachyonic matter [14], Phantom [15,16] and Chap- 
lygin gas [17] provide different time dependent functions 
for the equation of state [16,18-24]. There are also phe- 
nomenological models, parameterize the equation of state 
of dark energy in terms of redshift [25-27]. Here we ex- 
amine a simple phenomenological parameterization for 
the variable dark energy, proposed by Wetterich (2004) 
[28]. In this parameterization the variable dark energy 
is expressed in terms of three parameters. The first two 
parameters are the dark energy density, Q\ and dark en- 
ergy equation of state wq at the present time. The third 
parameter b is the bending parameter which can be ex- 
pressed in terms of the fraction of dark energy at early 
universe. The equation of state of this model depends on 
b and wq as: 
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Using the continuity equation, the density of dark energy 
changes with the rcdshift as: 



(2) 



where w(z;b,wo) = Wo/ [1 + 61n(l + z)\. In this paper 
we compare this model in a flat universe with the cosmo- 
logical observations such as SNIa, CMB shift parameter 
and Large Scale Structure (LSS) and constrain the pa- 
rameters of the model. 

The organization of the paper is as follows: In Sec. II 
we study the effect of the parameters of the dark energy 
model on the age of Universe, comoving distance, comov- 
ing volume element and the variation of angular size by 
the redshift [29]. In Sec. Ill we use the Gold sample of 
Supernova Type la data [30] to constrain the parameters 
of the model. In Sec. IV the position of the observed 
acoustic angular scale on CMB is compared with that of 
quintessence model. In Sec. V with the combined SNIa+ 
CMB+ SDSS data we put better constraints on the pa- 
rameters of this model. Finally the age of the universe 
determined by this model is compared with the age of the 
old stars and old high rcdshift objects as the consistency 
test. The conclusions are given in Sec. VI. 



II. THE EFFECT OF VARIABLE DARK ENERGY 
ON THE GEOMETRICAL PARAMETERS 

In this section we study the geometrical effect of the 
dark energy on the observable parameters of the universe 
as itemized as follows: 

• comoving distance: The radial comoving distance 
of an object located at a given redshift z is one of 
the basis parameters of the cosmology. Using the 
null geodesies in the FRW metric, the comoving 
distance can be obtained by: 



r(z;b, w ) 



dz' 



o H(z';b,w )' 



(3) 



where H (z; &, wq) is the Hubble parameter and for 
the redshifts smaller than the radiation dominant 
epoch (z < z eq ) 1 it can be expressed in terms of 
Hubble parameter at the present time, Ho, the mat- 
ter and dark energy content of the universe as: 



H 2 (z; b, w ) = Hl[tt m (l + zf + fi A (l + z) 



3[1-H5(z)]l 



(4) 



By numerical integration of equation (3), the co- 
moving distance as a function of redshift can be 
obtained. The dependence of comoving distance 
in terms of redshift for various bending factors is 
shown in Figure 1. This diagram shows that the 
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FIG. 1. Comoving distance, r(z;b, wo) (in unit of c/Hq) as 
a function of redshift for various values of bending parameters. 

comoving distance is more sensitive to b at higher 
redshifts than lower redshifts. Increasing the bend- 
ing parameter leads the dark energy dominance at 
the higher redshifts which results a slow growth of 
r(z; b, wo). For redshifts around z = we can ex- 
pand equation (3) as: 



z - -z 2 (i + n x w ) + 



(5) 



One of the main applications of the comoving dis- 
tance is in the luminosity distance calculation. In 
the next section we use SNIa data as the standard 
candle in the cosmological scales to confine the pa- 
rameters of the dark energy model. 

• angular size: The apparent angular size of objects 
at the cosmological scales is another observable 
that can be affected by the dark energy model. If 
D is the physical size of an object that subtends an 
angle 9 to the observer, for small 9 we have: 



D = d A 9 



(6) 



where d A = r(z; b, Wo)/(l + z) is called the angular 
diameter distance. One of the main applications of 
equation (6) is on measurement of matter content 
of the universe by observing the apparent angu- 
lar size of acoustic peak on CMB map and bary- 
onic acoustic peak at lower redshifts. A variable 
dark energy can change the comoving distance to 
the observer and consequently the apparent size of 
the acoustic peak. So measurement of the angular 
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FIG. 2. Alcock-Paczynski test, compares Az/A# normal- 
ized to the case of ACDM model as a function of redshift. 

size of objects in various redshifts (so-called Alcock- 
Paczynski test) can probe the variable dark energy 
[29]. The variation of apparent angular size A9 in 
terms variations of redshift Az can be determined 
as: 

Az _ H(z;b, w )r(z;b, w ) 

A6~ U 

Figure 2 shows Az/ AO in terms of redshift, nor- 
malized to the case with b = (ACDM model). 
According to Figure 2, the variation of the appar- 
ent angular size is sensitive to the bending param- 
eter at higher redshifts. The advantage of Alcock- 
Paczynski test is that we need standard ruler in 
the universe instead of the standard candle. Re- 
cently a possible Alcock-Paczynski test for mea- 
suring the dark energy parameters has been pro- 
posed by high-z observation of the power spectrum 
of large scale structures at 21cm wavelength. Ob- 
servation of 21cm fluctuations will enable us to de- 
termine the angular diameter and the Hubble con- 
stant [31,32]. The Ly-a forest of close QSO pairs 
may also measure £l\ and its variation with time 
[33]. 

• comoving volume element: The other geometrical 
parameter is the comoving volume element which 
is the basis of number-count tests, such as lensed 
quasars, galaxies, or clusters of galaxies. The co- 
moving volume element in terms of comoving dis- 
tance and Hubble parameters is given by: 

f(z;b,w ) = -^-r = r 2 {z- : b,w )/H(z;b 7 w ). (8) 
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FIG. 3. The comoving volume element in terms of redshift 
for various bending parameters. Increasing the bending pa- 
rameters makes the position of maximum value of volume 
element shifts to the smaller redshifts. 

As shown in Figure 3, the comoving volume clement 
reaches to its maximum value around z ~ 2 and for 
larger bending parameters the position of the peak 
shifts to the smaller redshifts. 

• age of the universe: The age of universe is another 
observable parameter that can be used to constrain 
the parameters of dark energy models. Studies on 
the old stars [34] suggests an age of 13!^ Gyr for 
the universe. Richer et. al. [35] and Hansen et. al. 
[36] also proposed an age of 12.7 ± 0.7 Gyr, by us- 
ing the white dwarf cooling sequence method. For a 
detailed review on the cosmic age see [5] . The uni- 
verse age in varying dark energy models depends 
on the parameters of the model. Integrating the 
age of universe from the beginning of universe up 
to now, we obtain the age of universe as follows: 

ft° r°° dz 

t (b,w )= / dt= -, (9) 

Jo Jo (1 + z)H{z; 6, wo) 

Similarly to the effect of bending parameter on the 
comoving distance, in the case of age of the uni- 
verse, increasing the bending parameter makes a 
shorter age for the universe. Figure 4 shows the 
effect of bending parameter on the age of universe. 
Here we show the variation of H Q t as a function 
of b for a typical values of cosmological parameters 
(e.g. h = 0.65, n m = 0.27 and w Q = -1.0). 
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FIG. 4. Age times Hubble constant as the present time 
(Hoto), as a function of b in a flat universe with the parameters 
of n m = 0.3, h — 0.65 and wo = —1.0. Increasing bending 
parameters makes a shorter age for the universe. 
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FIG. 5. Comparison of the distance modulus of the SNIa 
Gold sample in terms of redshift with the dark energy model. 
Solid line shows the best fit (Xmin/Nd.o.f = 1.13) with the cor- 
responding parameters of wo = — I.90I3 29, = O.Olj^oi 
and b = 6.00tl'fm in Iff level of confidence. 



III. TEST BY SUPERNOVA TYPE IA: GOLD 
SAMPLE 
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(12) 



The Supernova Type la experiments provided the main 
evidence for the existence of dark energy in the frame- 
work of standard cosmology. Since 1995 two teams of the 
High-Z Supernova Search and the Supernova Cosmology 
Project have been discovered several type la supernova 
candidates at the high redshifts [20,37]. Recently Riess et 
al. (2004) announced the discovery of 16 type la super- 
nova with the Hubble Space Telescope. This new sample 
includes 6 of the 7 most distant (z > 1.25) type la super- 
novas. They determined the luminosity distance of these 
supernovas and with the previously reported algorithms, 
obtained a uniform Gold sample of type la supernovas, 
containing 157 objects [30,38,39]. In this section we com- 
pare the distance modulus of the Gold sample data with 
that theoretically derived from the dark energy model. 
The distance modulus (// = m — M) in terms of redshift 
and parameters of model is given by: 



m — M = 5 log Dl (z; b, wq) + 25, 



(10) 



where M is the absolute magnitude, Dl is the luminosity 
distance in Mpc and m is the corrected apparent mag- 
nitude, including reddening, K correction etc. For a flat 
and homogeneous cosmological model the luminosity dis- 
tance can be obtained by: 



D L (z;b,w ) = (1 + z) 



dz' 



o H(z';b,w ) 



(11) 



The comparison between the observed and theoretical 
distance modulus is given by x 2 as follows: 



where o~i is the error bar of the observed distance modulus 
for each Supernova candidate. The best fit values for the 
model parameters are wq = — 1. 90^29' = 0.01+g'oi 
and b = 6.00±£j$ with X 2 nin /N d . . } = 1.13 at la level 
of confidence. Figure 5 compares the distance modulus 
of the observed SNIa Gold sample in terms of redshift 
and the best fit from the dark energy model. We see 
clearly that the fit values in this model are evidently 
different from those of ACDM (the WMAP results for 



ACDM models are: h = 0.71 



+0.04 

-0.03 



and fL 



0.27 



+0.04 
-0.04 



[5,40]). We also compare our result with that of Riess et 
al. (2004), putting b = we should recover their result. 
Figure 6 indicates the confidence contours of la, 2a and 
3cr in the (f2 m , wq) plane for the case of b = 0, marginal- 
ized over h, shows a good agreement with that of Riess 
et al. (2004). 

By substituting the cosmological parameters derived 
from the SNIa fit in equation (9), we obtain the age of 
universe about 13.45 Gry, which is in good agreement 
with the age of old stars. 



IV. THE LOCATION OF ACOUSTIC PEAK ON 
CMB MAP 

In this section we use the CMB data from the WMAP 
experiment to put additional constraints on the parame- 
ters of the dark energy model [41]. The statistical prop- 
erties of the temperature fluctuations on CMB is given 
by a two point correlation function. In the isotropic uni- 
verse, the two point correlation function C(j) depends 
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FIG. 6. Joint confidence intervals for Q m and wo for the 
case of b — with la (solid-line), 2a (dashed- line) and 3a 
(long dashed-line) confidence level. This result is in good 
agreement with that of Riess et al. (2004). 
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FIG. 7. Dependence of acoustic angular scale l A on the 
bending parameter, b and Wq for the three cases of l A — 290 
(solid-line), 300 (dashed-line) and 320 (dashed-dotted line). 



only on the angle between the two vectors (7) connect- 
ing the observer to the last scattering surface and it can 
be expanded into a Legendre polynomials as: 



C{l) 



1 OO 

— Y(2l 



l)QPz(cos 7 ), 



(13) 



1=2 



where C/ is the widely used angular power spectrum. The 
summation over I starts from 2 because the / = term is 
the monopole which in the case of statistical isotropy the 
monopole is constant, and it can be subtracted out. The 
dipole I = 1 is due to the local motion of the observer 
with respect to the last scattering surface and can be 
subtracted out as well. 

The relevant parameter in the spectrum of CMB which 
can determine the geometry and matter content of uni- 
verse is the position of the apparent acoustic peak. The 
acoustic peak corresponds to the Jeans length of photon- 
baryon structures at the last scattering surface some 
~ 379 Kyr after the Big Bang [5]. The position of the 
acoustic peak in Legendre polynomial space relates to its 
apparent angular size, 9 a through I a = it /6 a- The ap- 
parent angular size of acoustic peak in a flat universe can 
be obtained by dividing of the comoving sound horizon at 
the decoupling epoch r s (zdec) to the comoving distance 
of observer to the last scattering surface r(zdec) as: 



Ts (.Zdec) 
r{z d ec) ' 



(14) 



where we take the redshift of the decoupling at Zd ec = 
1089 [42]. The sound horizon corresponds to a distance 



that a perturbation of pressure can travel from the be- 
ginning of universe up to the decoupling area. The size 
of sound horizon at numerator of equation (14) can be 
obtained by: 



r s (z d e C ;b, w ) 



v s (z) 



H(z;b,w ) 



-dz 



(15) 



where v s (z)~ 2 = 3 + 9/4 x pb(z) / ' p r {z) is the sound ve- 
locity in the unit of speed of light from the big bang up 
to the last scattering surface [22,42]. The denominator 
of equation (14), r(zd ec ), comoving distance to the last 
scattering surface is also given by equation (3). 

Changing parameters of the dark energy can shift the 
size of apparent acoustic peak and the position of Ia in 
the power spectrum. Here we plot the dependence of l A 
on b and wq for a typical values of cosmological param- 
eters (Figure 7). It is seen that increasing b makes a 
shorter comoving distance to the last scattering surface 
(see Figure 1) and subsequently results in a larger acous- 
tic size or smaller I a- By a similar argument as in the 
case of b, shifting wq towards zero makes smaller I a- 

In order to compare the angular size of acoustic peak 
from this model with the observation we use the shift 
parameter as [43]: 



R 



dz 



E(z]b,w ) 



(16) 



where E(z\ b, wq) = H(z\ b, wq)/Hq. The shift parameter 
is proportional to the size of acoustic peak to that of fiat 



pure-CDM, A = model, (R oc 



fiat 



l fl A at /l A )- The 



observational result of CMB experiments correspond a 
shift parameter of R = 1.716 ± 0.062 (given by WMAP, 
CBI, ACBAR) [5,40]. One of the advantages of using 
the parameter R is that it is independent of the Hubble 
constant. The best ht values for the dark energy model 
using the combined CMB and SNIa observations, mini- 
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results in: f2 T - 



0.37 



-2.5±°f 5 With Xmin/Nd.o.f 



+0.08 
-0.07J 

1.11 



with (ler) level of confidence (see Table I). The age of 
universe from these parameters is 13.47 Gyr. 

In the next section we use the additional LSS data from 
SDSS and combine them with that of CMB and SNIa, to 
put more constraint on the parameters of the dark energy 
model. 



V. LSS COMBINED ANALYSIS WITH THE CMB 
AND SNIA 

Recent observations of large scale correlation function 
measured from the spectroscopic sample of 46,748 Lu- 
minous Red Galaxies (LRG) by the Sloan Digital Sky 
Survey shows a well detected peak around 100 Mpc hr 1 . 

This peak is an excellent match with the pre- 
dicted shape and the location of the imprint of the 
recombination-epoch acoustic oscillation on the low- 
redshift clustering matter [44]. For a flat universe we 
can construct parameter A as follows: 



A = y/UmE(zi;b,w ) 



-1/3 



dz 



1 



zi Jo E(z;b, w ) 



2/3 



(17) 



where E(z; b, wo) = H(z; b, wq)/Hq. Here A is also inde- 
pendent of Hq- We use the robust constraint on the dark 
energy model using the value of A = 0.469 ± 0.017 from 
the LRG observation at Z\ = 0.35 [44]. 

In what follows we perform a combined analysis of 
SNIa, CMB and SDSS to constrain the parameters of 
dark energy model. The combined x 2 from each experi- 
ment is as follows: 



X 



XSNIa 



2 i 2 

XCMB + XSDSSi 



(18) 



where XsNia * s obtained from equation (12), comparing 
the distance modulus of the SNIa candidates from the 
Gold sample with that of theoretical dark energy model. 
Xcmb ^ s a ^ so calculated by comparison of the observed 
shift parameter by WMAP with the model through equa- 
tion (16) and finally Xsdss * s obtained from equation 
(17). The marginalized likelihood functions (£oc e~ x / 2 ) 
[45] in terms of the parameters of the model for three 
cases of (i) fitting model with Supernova data, (ii) with 
SNIa+CMB data and (in) considering all three experi- 
ments of CMB+SNIa+SDSS are shown in Figure 8. 



TABLE I. The best values for the parameters of variable 
dark energy as Q m , b, wo with the corresponding age for 
the universe from the fitting with the SNIa, SNIa+CMB and 
SNIa+CMB+SDSS experiments with la and 2cr confidence 
level. 



Observation 
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age 
(Gyr) 


SNIa 


o oi +0 - 51 

u - ul -o.oi 

o oi +0 - 56 

u - ul -o.oi 


fi nn +7 - 35 

°- uu -6.00 
a nn+17.42 


-1 qn +0 - 75 

l.»U_ 3 2 g 

-I.90i7.2a 


13.45 


SNIa+CMB 


n 07+0. 08 
u -°' -0.07 

S7+ - 18 

u -°' -0.15 


7H+5.25 
z - ' u -1.65 

9 7n+11.85 
z - ' u -2.70 


-2.50t°; 9 8 » 
-2.50±i;S° 


13.47 


SNIa+CMB 
+SDSS 


n 97+0.02 
u - z '-0.02 

n 27+ 04 
u ' z ' -0.03 


1 or+1.65 
1 -°°-0.90 

1.35«; 3 3 ° 


-l-45t :SS 
-1.451°;" 


14.09 



O — 97+ 02 

J L m — u.z, 1 _Q.02' 

1.45tQ5Q at ler confidence 



The best fit values for the model parameters by 
marginalizing on the remained ones are 

b = 1.35±o;|a and w a = 
level with Xmin/Nd.o.f = 1.11. Table I indicates the cor- 
responding values for the cosmological parameters from 
this fitting with one and two a level of confidence. The 
joint confidence contours in the (wQ,£l m ), (wo,b) and 
(&, fl m ) planes are shown in Figures 9, 10 and 11. Com- 
paring Figure 6 which corresponds to the fitting with 
SNIa, prior b = with the solid contour in Figure 9, 
shows that adding the b parameter increases the degen- 
eracy between O m and Wo- 

We repeat the same analysis for the special case, fix- 
ing wo = —1.0. The best fit values for the parameters 
of the model in this case obtain as f2 m = 0.28^0 02) 
b = 0.24t°;23 at la confidence level with x min / N d.o.f = 
1.12. For the case of using the result of HST-Key project 
(Hubble parameter h = 0.71 ± 0.07), the marginalized 
likelihood function in this special case, where we have 
two free parameters of O m and b are shown in Figure 
12. The best fit values for the model parameters are: 
Q m = 0.20±g;{£, b = O.OOtom at lcr confidence level 
with Xmin/Nd.o.f = 1.75. In this case we have almost 
ACDM model with no variation in the equation of state 
of dark energy. 

Finally we do the consistency test, comparing the age 
of universe derived from this model with the age of old 
stars and the age of Old High Redshift Galaxies (OHRG) 
in various redshifts. One of the reasons for the exis- 
tence of the dark energy is the problem of "age cri- 
sis", i.e. the universe without cosmological constant is 



younger than its constituents [46]. From Table I, we see 
that the age of universe from the combined analysis of 
SNIa+CMB+SDSS is 14.09 Gyr which is in agreement 
with the age of old stars [34] . Here we take three OHRG 
for comparison with the dark energy model, namely the 
LBDS 53W091, a 3.5-Gyr-old radio galaxy at z = 1.55 
[46], the LBDS 53W069, a 4.0-Gyr-old radio galaxy at 
z = 1.43 [47] and a quasar, APM 08279+5255 at z = 3.91 
with an age of t = 2.1+ ) '5 ) Gyr [48]. The later one has once 
again led to the " age crisis" . An interesting point about 
this quasar is that it cannot be accommodated in the 
ACDM model [49]. To quantify the age consistency test 
we introduce the expression r as: 



t(z;b,w Q ) t(z;b,w )H 



^obs 



tobsHo 



(19) 



where t(z) is the age of universe which can be obtained 
from the equation (9) and t b s is an estimation for the age 
of old cosmological object. In order to have a compatible 
age for the objects, we should have r > 1. 

Table II shows the value of t for three mentioned 
OHRG. Various observational constraints on the param- 
eters of the dark energy model from SNIa, CMB, LSS 
and their combinations, results an age for the universe 
more than the age of LBDS 53W069 and LBDS 53W091, 
while APM 08279 + 5255 at z = 3.91 is older than the 
age of universe. Only in the case that we fix wq = —1 
and use the SNIa+HST constraints, we obtain r = 1.26, 
a compatible age for the universe with that of Quasar. 

VI. CONCLUSION 
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In this work we examined a parameterized quintessence 
model proposed by Wetterich (2004) [28] in a flat uni- 
verse with the low, intermediate and high redshift obser- 
vations. The effect of this model on the age of universe, 
radial comoving distance, comoving volume clement and 
the variation of apparent size of objects with the redshift 
(Alcock-Paczynski test) have been studied. 

In order to constrain the parameters of model we 
used the Gold sample SNIa data. The supernova anal- 
ysis results in a large degeneracy between the param- 
eters of the model. To improve the analysis, we com- 
bined analysis of SNIa and CMB shift parameter. Fi- 
nally we used the results of recently observed baryonic 
peak at the Large Scale Structure, combined with the 
two former experiments to confine the parameters of 
the dark energy model. The fitting parameters from 
the joint analysis of SNIa+CMB + SDSS marginaliz 
ing on the remained ones results in: fi ro 



n 97 +0.02 
u - z ' -0.02' 



b = 1.35io;go and w o = -l- 45 -o!60 at 1(7 confidence level 
with Xmin/Nd.o.f = Lll. The best fit value for the equa- 
tion of state of the dark energy leads to wo < — 1 which 
violates the strong energy condition. From the quan- 
tum field theory point of view, exotic models like scalar 




FIG. 8. Marginalized likelihood functions of three cosmo- 
logical parameters. The solid line is the likelihood function 
fitting the model with SNIa data, the long dashed-line with 
the joint SNIa+CMB and the dashed-line corresponds to the 
fitting with SNIa+CMB+SDSS data. The intersections of 
the curves with the horizontal solid and dashed lines give the 
bounds with la and 2a level of confidence respectively. 



TABLE II. The value of r for three high redshift objects, 
using the parameters of the dark energy model from the best 
fit. 



Observation 


LBDS 53W069 
z = 1.43 


LBDS 53W091 
2 = 1.55 


APM 
08279 + 5255 
z = 3.91 


SNIa 
, (w = -1) 


1.13 


1.21 


0.78 


SNIa+HST 

, (w = -1) 


1.67 


1.81 


1.26 


SNIa+CMB 
+SDSS 

, (w = -1) 


1.18 


1.27 


0.82 


SNIa+CMB 
+SDSS 
+HST, 
(wo = -1) 


1.29 


1.38 


0.89 


SNIa 


1.00 


1.05 


0.65 


SNIa+CMB 
+SDSS 


1.09 


1.17 


0.75 
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FIG. 10. Joint confidence intervals for b and wo from fit- 
ting with the SNIa (solid line) , SNIa+CMB (dashed-line) and 
SNIa+CMB+SDSS (long dashed-line) data with la level of 
confidence. 
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FIG. 9. Joint confidence intervals for fi m and w from fit- ting with the SNIa (solid line), SNIa+CMB (dashed-line) and 

ting with the SNIa (solid line), SNIa+CMB (dashed-line) and SNIa+CMB+SDSS (long dashed-line) data with la level of 

SNIa+CMB+SDSS (long dashed-line) data with la level of confidence, 
confidence. 




b 



FIG. 12. Marginalized likelihood functions for two cosmo- 
logical parameters of Q m and b in the special case of fixing 
wo = -1 and using the H = 71.0 ± 7.0 from the HST-Key 
project. The solid line corresponds to the Likelihood function 
of the parameters for the case of fitting dark energy model 
with the SNIa data, the dashed line corresponds to the com- 
bined SNIa+ CMB data and dotted-dashed line is for com- 
bined SNIa+CMB+SDSS data. The intersections with the 
horizontal solid and dashed lines give the bounds for la and 
2cr confidence level respectively. 

field with negative kinetic energy can provide w < — 1 
[50,51] (theoretical attempts for w < —1 can be found in 
[52-56]) 

We also did the age test, comparing the age of old 
stars and old high redshift galaxies with the age that we 
obtained based on the dark energy model. The age of 
the universe from the best fit parameters of the model, 
results in an age of 14.09 Gyr for the universe which is 
in agreement with the age of old stars. We also chose 
two high redshift radio galaxies at z = 1.55 and z = 
1.43 with a quasar at z = 3.91. The two first objects 
were consistent with the age of the universe, meaning 
that there were younger than the age of universe at the 
corresponding redshifts while the latter one was older 
than the age of universe. The age of APM 08279 + 5255 



quasar as the "age crisis" was not compatible with the 
age of universe in this dark energy model. Only in the 
case that we fixed Wo = —1 and took Hq = 71.0 ± 7.0 a 
compatible age for the universe with that of quasar has 
been obtained. 
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